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Abstract — A Physical layer Network Coding (PNC) scheme is 
proposed for the J-f -user wireless Multiple Access Relay Channel 
(MARC), in which K source nodes transmit their messages to 
the destination node D with the help of a relay node R. The 
proposed PNC scheme involves two transmission phases: (i) Phase 
1 during which the source nodes transmit, the relay node and 
the destination node receive and (ii) Phase 2 during which the 
source nodes and the relay node transmit, and the destination 
node receives. At the end of Phase 1, the relay node decodes 
the messages of the source nodes and during Phase 2 transmits a 
many-to-one function of the decoded messages. Wireless networks 
in which the relay node decodes, suffer from loss of diversity 
order if the decoder at the destination is not chosen properly. 
A novel decoder is proposed for the PNC scheme, which offers 
the maximum possible diversity order of 2, for a proper choice 
of certain parameters and the network coding map. Specifically, 
the network coding map used at the relay is chosen to be a K- 
dimensional Latin Hypercube, in order to ensure the maximum 
diversity order of 2. Also, it is shown that the proposed decoder 
can be implemented by a fast decoding algorithm. Simulation 
results presented for the 3-user MARC show that the proposed 
scheme offers a large gain over the existing scheme for the if -user 
MARC. 

I. Background and Preliminaries 

We consider the K-usei Multiple Access Relay Channel 
(MARC) shown in Fig. [T] Source nodes 5*1 , 52 , . . . , Sk want 
to transmit messages to the destination node D with the help 
of the relay node R. All the nodes are assumed to have half- 
duplex constraint, i.e., the nodes cannot transmit and receive 
simultaneously in the same frequency band. In a K-nser 
MARC, the maximum diversity order obtainable is two, since 
in addition to the presence of direct links, communication 
paths exist from the source nodes to the destination node D 
through the relay node R. 

A. Background 

In a fading scenario. Multiple Input Multiple Output 
(MIMO) antenna systems provide gain in terms of spatial 
diversity. However, in many practical scenarios, it is difficult 
to place multiple collocated antennas in a single terminal. An 
attractive alternative to obtain diversity gain without using 
multiple antennas, is the utilization of intermediate relay 
nodes which aid the transmission from the source nodes to 
the destination nodes. In order to exploit the presence of 
intermediate relay nodes to obtain diversity gain, the source 
nodes need to convey their messages to the relay nodes. Due 
to the superposition nature of the wireless channel, if the 
source nodes transmit simultaneously in the same frequency 
band, interference occurs at the relay nodes. A loss of spectral 



efficiency results, if the source nodes transmit in orthogonal 
time/frequency slots. A solution to this problem is the use of 
physical layer network coding, first introduced in |[T|, in which 
the nodes are allowed to transmit simultaneously. Physical 
layer Network Coding (PNC) has been shown to outperform 
traditional schemes which involve orthogonal transmissions 
|T|- iTl. So far, most of the works on PNC have mainly 
focussed only on the two-way relay channel. In our recent 
work Q, we proposed a scheme based on PNC for the two 
user MARC. In this paper, we present the generalization of 
the scheme proposed in |^ for the K-nsei MARC. 

For the MARC, a Complex Field Network Coding (CFNC) 
scheme was proposed in fT|. The CFNC scheme, like the PNC 
scheme, avoids the loss of spectral efficiency, by making the 
source nodes transmit simultaneously. But the major difference 
between the CFNC scheme and the proposed PNC scheme is 
that during the relaying phase, the CFNC scheme uses a signal 
set of size at R, whereas the proposed PNC scheme uses 
a signal set of size M, where M is the size of the signal set 
used at the source nodes. 

As noted in |[7), if the relay node transmits a many -to- 
one function of the estimates of the messages transmitted by 
the source nodes and minimum squared Euclidean distance 
decoder is employed at D, a loss of diversity order results. This 
problem of loss of diversity order due to error propagation, 
is encountered in many other wireless scenarios as well and 
various solutions have been proposed to avoid this problem. 
Cyclic Redundancy Check bits are used so that the nodes 
forward only those packets which are decoded correctly |8|. 
Some works assume the knowledge of all the instantaneous 
fade coefficients or error probabilities associated with the 
intermediate nodes at the destination node, with the decoder 
at the destination using this knowledge to ensure full diversity 
Q, fiQ\ - The CFNC scheme proposed in fl\, uses a scaling 
factor at the relay node which depends on the instantaneous 
fade coefficients associated with the hnks from the source 
nodes to the relay node, with the scaling factor indicated to the 
destination using pilot symbols. The proposed scheme does not 
suffer from the disadvantages of any of the above methods, yet 
ensures the maximum possible diversity order. This is achieved 
by means of an efficient choice of the transmission scheme and 
the use of a novel decoder at the destination D. 

For the proposed PNC scheme, making the source nodes 
also transmit during the relaying phase, combined with a novel 
decoder ensures the maximum possible diversity order of two. 
Furthermore, if certain parameters are chosen properly, the 



proposed decoder for the PNC scheme can be implemented 
with a decoding complexity order same as that of the CFNC 
scheme proposed in |7|. 

For the two-way relay channel, the network coding maps 
used at the relay node need to form a mathematical structure 
called Latin Squares, for ensuring unique decodability at the 
end nodes [11 1. The structural properties of Latin Squares 
have been used to obtain the network coding maps in a 
two-way relay channel [llj- | ,13J . Interestingly, choosing the 
network coding map used at R to be a iC-dimensional Latin 
Hypercube, which is the generalization of the Latin Square to 
K dimensions, helps towards achieving the maximum diversity 
order of two for the K-usei MARC. 




(a) Phase 1 (b) Phase 2 

Fig. 1: The two user Multiple Access Relay Channel 
The main advantages of the proposed scheme over the 
CFNC scheme proposed in \JJ are summarized below: 

• In the CFNC scheme, R transmits a complex linear 
combination of the estimate of the messages transmitted 
by the source node and the signal set used at R during 
the relaying phase has points, where M is the 
size of the signal set used at the source nodes. The 
minimum distance of the signal set used during the 
relaying phase vanishes as K increases. In contrast, since 
the proposed PNC scheme uses a many-to-one map, the 
signal set used during the relaying phase has only M 
points. The minimum distance of the signal set used 
at R is more than that of the CFNC scheme and it 
remains the same irrespective of the number of source 
nodes K. Hence the proposed scheme performs better 
than the CFNC scheme. Simulation results presented for 
the 3-user MARC confirm that the proposed PNC scheme 
provides a large gain over the CFNC scheme. 

• In the CFNC scheme, R uses a scaling factor which is a 
function of the instantaneous fade coefficients associated 
with the links from the source nodes to the relay node, 
which needs to be indicated to D using pilot symbols. 
Since the proposed PNC scheme does not involve any 
such scaling factor, there is no need of such pilot symbols. 

Notations: Throughout, vectors are denoted by bold lower case 
letters and matrices are denoted by bold capital letters. The 
set of complex numbers is denoted by C. CAf{0, cr^) denotes 
a circularly symmetric complex Gaussian random variable 
with mean zero and variance cr^ and J\f{0,a'^) denotes a 



real Gaussian random variable with mean zero and variance 
a^. For a matrix A, A'^ and A* denotes its transpose and 
conjugate transpose respectively. For a complex number x, 
X* denotes its conjugate and |a;| denotes its absolute value. 
For a vector v, ||v|| denotes its Euclidean norm. The total 
transmission energy at a node is assumed to be equal to Eg 
and all the additive noises are assumed to have a variance 
equal to 1. By SNR, we denote the transmission energy Eg. 
For a signal set S, AS denotes the difference signal set of 
S, AS = {x — x'\x,x' e S}. The all zero matrix of size 
n X n is denoted by 0„. The natural logarithm of x is denoted 
by log(x). E{X) denotes the expectation of X. Q[.] denotes 
the complementary CDF of the standard Gaussian random 
variable. 

B. Signal Model 

Throughout, a quasi-static fading scenario is assumed with 
the channel state information available only at the receivers. 
The source nodes want to transmit a binary vector of length A 
to the destination node. At each one of the source nodes, the 
binary vector is mapped onto a point from a M = 2^ point 
signal set denoted by S. Let : F2 — > 5 denote the mapping 
from bits to complex symbols used at the source nodes. 

The proposed PNC scheme involves two transmission 
phases: Phase 1 during which the source nodes simultaneously 
transmit and, R and D receive, followed by the Phase 2 during 
which the source nodes and R transmit to D. 

Phase 1: Let Xi = ^(s,) € 5, i € {1,2... ,K},s.i G 
denote the complex symbol the source node Si wants to 
convey to D. During Phase 1, the source node Si transmits a 
scaled version of Xi. The received signal at R and D during 
Phase 1 are respectively given by, 

K 

y_R = ^ hs^R^/E^sCiiXi + zr and 

i=l 
K 

ypi ^ hs^D\/E~saiXi + +ZD^ , (1) 

where ai E C,i E {1,2 ... , K} are constants and the additive 
noises zr and z^^ are assumed to be CAf{0, 1). All the fade 
coefficients are Rayleigh distributed, with the fade coefficient 
associated with the Si-R link hs-R ^ CA/'(0, cr|.^), and 
the fade coefficient associated with the Si-D link hs d ^ 
CN{Q,<jIo). 

Based on the received complex number yR, the relay 
node computes the Maximum Likelihood (ML) estimate of 
(xi, X2, . . . , xk) denoted by (xf , X2 , . ■ . , x^), i.e., 

K 

(if ,f J, . . . = arg ^ , niin _\yR-'^hs^R\/E^aix'i\. 

Phase 2: During Phase 2, the source node Si trans- 
mits a scaled version of Xi and R transmits xr — 
f{xf,X2,... ,x^), where / : S^ — 5 is a many-to-one 
function. The received signal at D during Phase 2 is given by. 



(2) 



where 6, g C are constants and the additive noise zjj^ is 
assumed to be CAf{0, 1). The fade coefficient associated with 
the R-D link is assumed to be CJ\f{0, aj^j-,). 

For the transmission energy at the source nodes to be equal 
to Eg, the constants and bi are chosen such that |aip + 

|6,f = i,v*e{i,2... 

From ([T]l and (|2|, the received complex numbers at D during 
the two phases can be written in vector form as, 



[vdi yn2] = \/E^ [liSin hs2D ■ ■ - hsj^D 'ird] 



a2X2 b2X2 



a-K-XK bKXK 
, Xr 



C(xi.x2,... ,!cK'^r) 
+ [zDj ZD2] ■ (3) 



The matrix C{xi,X2, ■ ■ ■ ,xk,xji) in (|3]l is referred 
to as the codeword matrix. The restriction of 
C{xi,X2, ■ ■ ■ ,xk,xr) to the first K rows, denoted by 

Crfxi, a;2, . . . , xk) is referred to as the restricted codeword 

r 1 r 

matrix. I.e., Cr(xi,x2,...,x^)= [,^1^1 62x2 ... fe^J ' 
The difference between any two restricted codeword matrices 
is referred as the restricted codeword difference matrix, i.e., 
the restricted codeword difference matrices are of the form 



Cr{Ax^,Ax2,... ,Axk) - ^^^^^^ ^^^^^ 

where Ax, e A5,Vi G {1,2, . . . ,K}. 

From the vector yn can also be written as. 



bxAxK 



1=1 



where is a {K + 1) x 2 matrix whose row is given by 
\ai hi] and all other entries are zeros. For the (if +1) x 2 matrix 
Wr, the {K + 1)* row is given by [0 1] and all other entries 
are zeros. The matrices Wi, i G {1,2... , K} and Wr are 
referred to as the weight matrices. 

The contributions and organization of the paper are as 
follows: A novel decoder for the proposed PNC scheme is 
presented in Section 11 A. In Section II B, it is shown that 
the decoder presented in Section II A achieves the maximum 
diversity order of two if the following two conditions are 
satisfied: (i) the map / used at the relay node forms a K- 
dimensional Latin Hypercube and (ii) the constants a-i and 
bi,i G {1,2... , K} are such that every 2x2 square submatrix 
of the restricted codeword difference matrices have rank two 
when Axi takes non-zero values. In Section III, the condition 
under which the proposed decoder admits fast decoding is 
obtained. It is shown that when at least one of the weight 
matrices W; is Hurwitz-Radon orthogonal with Wr, the 
proposed decoder admits fast decoding, with the decoding 
complexity order same as that of the CFNC scheme proposed 



in |[7|. Simulation results which show that the proposed 
PNC scheme provides large gain over the CFNC scheme are 
presented in Section IV. 

II. A Novel Decoder for the Proposed PNC Scheme 
AND ITS Diversity Analysis 

In the following subsection, a novel decoder for the pro- 
posed PNC scheme is presented. 

A. A Novel Decoder for the Proposed PNC Scheme 

When D uses the minimum squared Euclidean distance 
decoder given by, 

r ^ 

( , if , . . . , £^ ) = arg min < Iud, - hs, d \/E^ai Xi\'^ 



+ lfD2 ^ hs^D\/ Esbi Xi - Hrd \/ Esf{xi ,X2,. . . , Xk)\ 

a loss of diversity order results, since this decoder does not 
consider the possibility of decoding errors at the relay node. 
Alternatively, we propose a novel decoder given by. 



: are mm 



{mi (xi, X2, . . . , Xk) , 



\og{SNR)+m2{xi,X2,... ,xk)}, (4) 

where the metrics nii and ni2 are given in (|5]l and (|6]) 
respectively, at the top of the next page. 

The idea behind the choice of this decoder is as fol- 
lows: The optimal ML decoding metric at D is equal to 
mi{xi,X2,--- ,xk), when the relay transmits the correct 
network-coded symbol. The relay transmits a wrong network- 
coded symbol, independent of {xi,X2, ■ ■ ■ ,xk), if the joint 
ML estimate at the relay {x^,X2,--- is such that 

XR = /(if,i|^,... ,x^) ^ f{xi,X2,-.- ,xk)- Under this 
condition, the optimal ML decision metric at D is given 
by m2{xi,X2, ■ ■ ■ tX^)- At high SNR, the relay transmits a 
wrong network-coded symbol with a probability which is pro- 
portional to glij^ - Hence, to the metric TO2(a;i, 2^2, ■ ■ ■ ■, xk), 
we add a correction factor of \og{SNR) and the minimum of 
TOi(xi,X2, . . . :Xk) and \og{SNR) + m2{xi,X2, ■ ■ ■ ,xk) is 
taken to be the decoding metric at D. 

The CFNC scheme proposed in |7| uses the minimum 
squared Euclidean distance decoder, which has a decoding 
complexity of 0{M^). Since the decoder given in Q involves 
minimization over K+1 variables xi,X2, ■ ■ ■ xk and xr., it ap- 
pears as though the decoding complexity order is C'(A/^+^). 
In Section III, it is shown that by properly choosing the 
constants a^'s and fo^'s, the decoding complexity order can be 
reduced to 0{M^) which is the same as that of the CFNC 
scheme. 

B. Diversity Analysis of the Proposed Decoder 

The following theorem gives a sufficient condition under 
which the decoder given in (|4]) offers maximum diversity order 
two. 

Theorem 1: For the proposed PNC scheme, the decoder 
given in Q offers maximum diversity order two if the fol- 
lowing two conditions are satisfied: 

1) The map / satisfies the condition. 



nil {xi; X2, 



- hs-D\/E^bi Xi - hao \/E^f(x-i_,X2, . . . , xk) 



(5) 



m2 X2, .... xk) — 


K 

VDi — hs-D\/E~sai Xi 

i=i 


J 

+ mill < 




K 

— ^2 ^S^D\/E~sbi Xi — IlRD \/E~sXr 




(6) 


ms {xi,X2, . . . .Xk) ^ 


K 

Vdi - hs-n\/E^ai Xi 
i=i 


J 

+ mm < 




K 
i=l 


Eghi Xi — hjijj \/ EsXji 




(7) 



f(xi,X2,... ,Xi-l,Xi,Xi + l,... ,xk) 

^ f{xi,X2, . . . ,Xi-l,x[,Xi+l, 



,xk), (8) 



for ^ a;-, for all i e {1, 2, . . . , K}. 
2) All 2x2 submatrices of the restricted codeword differ- 
ence matrices Cr(Axi, Ax2, ■ ■ ■ , Axx) have rank two, 
VAa;i,Aa;2,... ,Axk ^ 0. 
It is easy to verify that a map / satisfying condition 1) 
above forms a Latin Hypercube of order AI and dimension 
K. 

Definition 1: | [T4) A Latin Hypercube of order M and 
dimension K is an array of dimension with the indices 
for the K dimensions as well as the entries filled in the 
array taking values from the symbol set {0, 1 . . . ,M— 1}. 
Every symbol occurs exactly once along each one of the K 
dimensions. 

The i* dimension of the Latin Hypercube represents the 
transmission Xi of the source node Si. For simplicity, the 
points of the M point signal set S are indexed by integers 
from to AI — 1. The entries filled in the Latin Hypercube 
represent the transmission of the relay node. 
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Fig. 2: An example of Latin Square of order 4 

A Latin Hypercube of dimension 2 is a Latin Square. Fig. 
|2] shows an example of a Latin Square of order 4. It can be 
seen from Fig. |2] that no two entries repeat in a row as well 
as a column. 



Cr(Aa;i, A2;2,Aa;3) = 



. It can be ver- 



Similarly, for a fixed value of row index and a fixed value of 
column index, the entries in the four pages are distinct. 

For the 3-user MARC, the following example gives a choice 
of a^'s and b/s for which Condition 2) given in Theorem 1 is 
satisfied. 

Example 1: For the 3-user MARC, choosing oi = 1, 

hi = 0, a2 = ^2 = ^, 03 = ^ and 63 - ^^^j, 
the restricted codeword difference matrices are of the form 
" Aa;i 

_^Ax3 --73^x-3j 
ified that when Azi, Aa;2 and AX3 take non-zero values, the 
rank of every 2x2 square submatrix of Ci.(Aa:i, Ax2, Aa;3) 
is two and hence condition 2) given in Theorem 1 is satisfied. 

III. A Fast Decoding Algorithm for the Proposed 
Decoder 

In this section, it is shown that if the constants a^'s and 
fej's are chosen properly, the decoder given in Q can be 
implemented using an efficient algorithm with a complexity 
order 0{AI^). Note that for the CFNC scheme proposed in 
I?], the decoding complexity order at D is 0{M^). 

Before the algorithm is presented, some notations are in- 
troduced. The points in the signal set S are denoted by 
Sj, 1 < i < M. 

From (|3]l, the vector can be written as, 

xi 

X2 







He 



:i 1 



Fig. 3: An example of Latin Cube of order 4 



The matrix Heq can be decomposed using QR decomposition 
as Heq = QR-, where Q is a 2 x 2 unitary matrix and R = 
[Ri R2] is a 2 x {K + 1) matrix, with Ri being upper- 
triangular of size 2x2 and R2 being a 2 x [K — 1) matrix. 
Let Tij denote the entry of R. 

Define yn = Q^yS = [VDi 1)02^- Also, let 



pl(Xi,X2, . . . 



A Latin Hypercube of dimension 3 is a Latin Cube. The 
three dimensions of a Latin Cube are referred to as rows, 
columns and pages. Fig. |3] shows an example of a Latin Cube 
of order 4. It can be verified from Fig. |3]that in each one of the 
four pages, no two entries repeat in a row as well as a column. 



i>2{x\,X2, ■ ■ ■ ,Xk) = ISdj ~ 'y^j2iXi\rEs 

i = 2 

- ''2(K + 1)/(^1, 2:2, . . . .,XK)xfW,\^ , 

K 

H{x2, ■ ■ ■ , Xk, xr) = |yr>2 - ^ T2iXi~/lfs - r2(K+l)a:ij\/iJri^- 



The following proposition gives a sufficient condition under 
which Algorithm 1 below implements the decoder given in 
©■ 

Algorithm 1 Decoding Algorithm used at D 



for 2-2 — 1 to M do 
for 13 = 1 to M do 



for ijc = 1 to M do 

X3 +- Si„ 



from Theorem 2, (15)). Hence, from ( [TO] i, it follows that, 

TO3(a;i, a;2, . . . , xk) = | j/Di - ^ riiXi\/E~s\'^ 

i=l 

+ min < \yD2 - 7-2ja;i\/£^ - f2(X+i)a;flV^P ^ • (H) 



Hence, 



a; A- ^ 
Find a 



Find : 



arg min {01 (xi , a;2 : 

: arg min {i^i{xi, X2, 
3=1 es 



■ , Xk) + <t>2{xi; X2, 

■ ■ ,xk)} 

, Xk,Xi{)} 



,Xk)} 



Find Xfi — arg min {03 (3:2 
if <l)i{xl,X2, ■ ■ ■ ..Xk) + 4>2(x\,X2, ■ ■ - Xk) < log{SNR) 

+ 4>l{xl,X2, . . . ,Xk) + 4'3(,X2, . . . ,Xk,Xr) 

then 

m{x2, ■ 

Xl(x2, . 

else 

"1(2:2, ■ ■ 



Xl{x2, 

end if 
end for 



, Xk) = <l>i(x\,X2, 
,xk) = xl 

, Xk) = <t>i(x\,X2, 

,xk) = x\ 



, XK) + <t>2{Xi,X2, . . . , Xk) 



,Xk)+<I>3{x2, . . . ,Xk,Xr) 

+ log{SNR) 



end for 
end for 

Find (£2, ■ • ■ ,xk) = arg 

(£f , X2 , ■ ■ ■ ; x'j^) = {xi(x2, 



mm 

,^k)^s 



K~l 
Xk):X2, 



m(x2 , . . . Xk) 
■ ■ ■ ,xk) 



Proposition 1: Algorithm implements the decoder in 
Q, if the constants a^'s and 6,;'s are such that the weight 
matrices Wi and Wr are Hurwitz-Radon (H-R) orthogonal, 
i.e., WiWr* + WrWi* = Ok+1. 

Proof: The decoding metric of the decoder given in (|4]) 
can be written as, 

min {mi (a;i, X2, . . . , xk) , log (SNR) + m2 (xi, X2^ ■ . ■ , xk)} 

(xi,X2,-.. ,xjf ) 

— min {mi (xi , X2., . . . , xk) ■. 

(xi,X2,... ,xj^) 

log (SNR) + mi (xi, X2, ■ ■ ■ , Xk) , log (SNR) + m2 {xi, X2, ■ ■ ■ , xk)}- 

— min {mi (xi , X2, . ■ . , xk) , (9) 

(xi,X2,... ,xjf ) 

log (SNR.) + m3 (xi, 2:2, ■ ■ • , Xk)}- 

where the metric 7713(2:1, X2, .. . ,xk) is given in (|7]l, at the 
top of the previous page. We have, 

7713(0:1, X2, .. . ,xk) = min {||yD^ - HeqxV^^lP} 

= min{||yD-Rx/B;||2}. (10) 

xrES 

Since Ri is upper triangular, r2i = 0. Also, the entry 
^i(A'+i) = 0, since Wa and Wr are H-R orthogonal (follows 



' A algorithm exactly similar to Algorithm 1 cati be used with the roles of 
Si and Si interchanged, if W; and Wh are H-R orthogonal. 



min 7773(0:1, 2:2, . . . ,xk) = min 01(2:1,2:2, . . . , xk) 
+ min (t)3ix2, ■ ■ ■ ,xk,xr). 

xrES 

From (|9|, the decoding metric can be written as, 

• ■ • ,xk) 

+ cj)2 (2:1,2:2, . . . ,2: if)} , 



mm < min {(/)i (2:1,2:2, ... ,2;^) 

(x2,...,xk)£S'<-^ \xies 



min cf)i (2:1,2:2, . . . ,xk) 

xi^S 



+ min fl!)3 (2:2, . . . ,xk,xr) + \og{SNR) } 
XRes J 



In Algorithm 1, inside the K—1 nested for loops, the values 
of 2:2, 2:3 . . . , Xk are fixed and the operations in lines 13, 14 
and 15 involve a complexity order 0{M). The operations from 
line 16 to line 23 involve constant complexity, independent of 
M. Hence the complexity order for executing the nested for 
loops from line 1 to line 29 is 0{M^). The operation in line 
30 involves a complexity order 0{AI^~^). Hence the overall 
complexity order of Algorithm 1 is 0{M^), which is the 
same as that of the CFNC scheme proposed in ||7|. 

Example 2: Continuing with Example 1, for the 3-user 
MARC, when ai = 1, 61 = 0, 02 = 4^, h 
and 63 — 



and 03 = 73, 
given by, Wi = 



the weight matrices are 



W3 



1 

V2 





n T 

and Wr 



W2 



75 
o' 

1 






It can 



be veiified that the matrices Wi and Wr are H-R orthogonal, 
i.e., WiWr* + WrWi* = O4. Hence, for this case. 
Algorithm 1 can be used to implement the decoder given 
in (|4]i. 

IV. Simulation Results 

Simulation results presented in this section compare the 
performance of the proposed PNC scheme with the CFNC 
scheme proposed in |l7T|, for the 3-user MARC. In all the 
simulation results presented, the values of the constants a/s 
and &i's are chosen to be the ones in Example 1 and 4-PSK 
signal set is used at the nodes. The Latin Cube given in Fig. 
[3] is used as the network coding map at the relay node. 

For the case when the variances of all the fading links are 
dB, the SNR Vs. Symbol Error Probability (SEP) plots are 
shown in Fig. |4] It can be seen from Fig. [4] that the PNC 
scheme performs better than the CFNC scheme and offers a 
large gain of 8 dB, when the SEP is 10^*. Fig. [s] shows a 



similar plot for the case when (t| — 10 dB and cr| — 
apjj = dB, where i G {1,2,3}. It can be seen from Fig. 
[5] that for this case, the PNC scheme offers a gain of nearly 
6 dB, when the SEP is 10^^. Fig. |6] shows the plots for the 
case when the R-D link is stronger than all other links, i.e, 
= ^ko = (^Id = dB,i e {1,2,3} and cr|^ = 10 
dB. For this case, the PNC scheme offers a large gain of about 
12 dB, when the SEP is lO^''. Also, it can be verified from 
the plots that the proposed decoder for the PNC scheme offers 
the maximum possible diversity order of two. 
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Fig. 4: SNR vs SEP plots for 4-PSK signal set for al.^ = 
c^liD = cTflD = dB. 



T 
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Fig. 5: SNR vs SEP plots for 4-PSK signal set for al^^ = 10 
dB, 0-5,0 = o-SD = o-flD = dB. 

V. Discussion 

A physical layer network coding scheme was proposed 
for the K-nsei Multiple Access Relay Channel. For the 
proposed scheme, a novel decoder was presented and it 
was shown that the decoder offers the maximum possible 
diversity order of two if the network coding map used at 
the relay forms a if-dimensional Latin Hypercube and every 
2x2 submatrix of a restricted codeword difference matrix 
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Fig. 6: SNR vs SEP plots for 4-PSK signal set for a^.j^ = 
o-|_o = dB, cr|o = 10 dB. 

Cr(Aa;i,a;2,.. . , Aa;^), Ax, 7^ 0,Vi G {1,2,... ,K}, has a 
rank two. Also, it was shown that the proposed decoder can 
be implemented using a fast decoding algorithm, if a weight 
matrix Wj is Hurwitz-Radon orthogonal with Wr, for some 
i E {1,2... , K}. The problem of finding the constants a/s 
and bi's which minimize the error probability in addition to 
ensuring maximum diversity order remains open. Extension 
of the proposed scheme for the case when there are multiple 
relay nodes is a possible direction for future work. 
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APPENDIX - Proof of Theorem 1 

Let H denote a particular realization of the fade coeffi- 
cients. Throughout the proof, the subscript in a probability 
expression indicates conditioning on the fade coefficients. For 
simplicity of notation, it is assumed that the variances of all 
the fading coefficients are one, but the result holds for other 
values as well. 

Let E denote an error event that the transmitted message 
if-tuple [xi. . . . , xk) is wrongly decoded at D. 

The probability of E conditioned on H 
in ( [T2I ), can be upper bounded as in ([T3j 
(IT2I1 and ( [T3| l are shown at the next 



Ph 
Ph 



R ^R 
^2 ' ■ 



= f{xi,X2,. 
7^ f{xi,X2,. 



given 
(eqns. 
page), 
and 

respectively 



denote the probabiUties that R transmits the correct and wrong 
network coded symbol during Phase 2, for a given H. Also, 
the probability Ph {{xi ,X2 , ■ . ■ ,Xk) ixi,X2, . ■ . ,xk) 
\ f(xi,X2,. ■ ■ ,x'^) = ,f {xi,X2,. ■ . ,xk)} and the 

probability Ph {{xf,x§, . . . ,x'^) / {x-i_,X2, ■ ■ ■ ,xk) 

\ f{x^,x^,... ,x'^)j^ f{x^,X2,... ,xk)} in ^ 

respectively denote the probabilities of E given that 
R transmitted the correct and wrong network coded 
symbol for a given H. Ph{E} can be upper bounded 
as in ([14]), where Pf/ |/(xf , if , . . . , ^) = 
denotes the probability that the network coded symbol 
transmitted by R is x'j^ ^ f{xi,X2,.-- ,xk) and the 



probability Ph { (xf , xf' , . . 



) 7^ {xi,X2, ■ . ■ ,xk) 



• i^k) — x'r\ the probability of E given 



that R transmits a;'^^, for a given H. Taking expectation of the 
terms in ([T4]i w.r.t we get ( fTSj ). 

The rest of the proof of Theorem 1 is presented in two parts 
as Lemma 1 and Lemma 2. In Lemma 1, it is shown that the 



probability 

I f{if,xl^: 

diversity order 
probability 



P { i^i , X2 , ■ ■ 



^^k) = fixi,X2, 



k) ixi,X2, ■ . ■ :Xk) 



two. Lemma 



,xk)} 
2 shows 



■ • J x^) — X [ 

~ i.R 



P{ix'^,X^,X^) ^ {XI,X2, 



has a 
that the 

,xk) 

c'fj} has a diversity order one. 
Since P {/(if , . . . ,x^)=x'j^} has a diversity order 
one. Lemma 1 and Lemma 2 together imply that P {E} has 
a diversity order two. 

Lemma 1: When the two conditions in the 
statement of Therorem 1 are satisfied, the 



D 

2 ' ■ 



■Xk) ^ {xi,X2, ■ . ■ ,xk) 
,xk)} has a diversity 



probability P | (xf , x 

I f{Xi ,X§,x'^) = f{xi,X2 

order two. 
Proof: 

Recall that the decoder used at D given in ([4]) in Section II 
A, involves computation of the metrics mi and TO2 defined 
in (|5]) and (|6]|. Under the condition that R transmitted the 
correct network coding symbol, a decoding error occurs at 
D only when mi{xi,X2, ■ ■ ■ ,xk) > mi{x[,X2, ■ ■ ■ ,x'k) 
or mi(xi,a;2, . . . ,xk) > \og{SNR) + 'm2{x\,x'2, ■ ■ ■ ^x'^) 
for some {x'i,X2, ■ ■ ■ ^x'^) ^ {xi,X2, ■ ■ ■ ,xk)- Hence, 
the probabihty P {{xi ,X2 , ■ ■ ■ ,x'j^) ^ {xi,x2, ■ ■ ■ ,xk) 
\ f{x^,x2,... ,x^) = f {xi,x2,... ,xk)} can be upper 
bounded as in ( [T6| ), which can be upper bounded using the 
union bound as in ( [T7| ) (eqns. ( [T6] l and ( [T7] i are given at the 
next page). 

The probability P {mi (xi, a;2, .. . ,xk) > 'm\{x'i,x'2, . . . ,x'j^) 
I f{xf, £2, . . . , xff) = / {xi,x2, ■ ■ ■ , Xk)} is equal to the Pair- 
wise Error Probability (PEP) of a space time coded 3x1 
collocated MISO system, with the codeword difference 
matrices of the space time code used at the transmitter 
being of the form ["i^^i ^'^^^^ ■ ■ "^'^"'^ a° 1^- 
where Axi — Xi — x^,i G {1,2,... ,K and 
AxR = f{xi,X2,-.- ,xk) - f{x[,X2,...,x'K)- When 
Axi 7^ 0, for at least two values of i e {1,2,... ,K}, 
these codeword difference matrices are of rank 2, 
other wise condition 1) given in the statement of 
Theorem 1) will be violated. When A.t^ 7^ and 
Axj = 0,Vi ^ jjhj G {1,2,... ,K}, the codeword 
difference matrices are full rank, otherwise condition 2) 
in the statement of Theorem 1 will be violated. Since 
the codeword difference matrices are full rank, the 
probability P{mi(a;i,a;2, . . . ,xk) > mi{x'i,X2, . . . ,x'jf) 

\ f(xi,X2,... ,xk) = f {xi,x2,... ,xk)} has a diversity 
order two ||T6). 

be a metric as defined in ( [T8] l, 
the next page. The probability 
. . . , Xk) > m2{x[,x'2, . . . , x'k) + log(5'iVP) 
. , Xk) = f {xi,x2, ■ ■ ■ , xk)} can be written in 
terms of the metrics nii and 1114 as in ( [T9] l, which can be 
upper bounded as in ( |20] i (eqns. ( [T9] l - ( |22] i are given in the 
next page). 



Let m4 
given in 

P{mi(a::i,X2, 



Let Axr = f{xi,X2, ■ ■ ■ ,xk) - x'j^, Ax^ 



x'. The 



probability P^ {mi(a;i, a;2, . . . ,XK)>rn4{x'i,X2,... ,x'k,x'ji) 
\ f{xi,x2, ■ ■ ■ ,xk) = f{xi,x2, ■ ■ ■ ,xk)} can be written in 
terms of the additive noise zj^^ and zu^, as given in ( |2T] i. 

Let xi ^ Yl,f=ihs^Dai\/EsAxi, X2 = {J2iLiif^SiDbt + 
h]i£iAxii)\/Ers- Also, let zd = [^Di zjj^] and x = 
[xi X2^ . Then ( pTj i can be simplified as in ( |22] l, 
where w ~ -\/2Pe{zQ^}, is distributed according to 
A/^(0, 1). In terms of the Q function, the probability 



Ph\w< 



- log(SNR) 



i} 



_ y^ii II — J 111 ( p2] i can be written as 

Q ' + ^ • Note that ||x|| depends on the fade 

coefficients. To complete the proof, it suffices to show that 



Ph{E} = Ph |(xf , ... , x^) 7^ (xi, X2, . . . ,xk) \ /(xf , x^, . . . , x^) = / (xi, X2, . . . , Xft:)| Pjj |/(xf , x^, . . . , x^) = / (xi, X2, . . . , Xft:)| 

+ f";/ {(xf.xf , . . . ,x^) 7^ (xi,X2, . . . , xk) I /(xf .xl*, . . . ,x^) 7^ / (xi,X2, . . . ,xk)} Ph {/(xf ,xf , . . . ,x^) 7^ / (xi,X2, . . . ,xk)^ 

(12) 

< Ph {(xf , a J, . . . , x^) / (xi,X2, . . . , Xftr) | /(xf ,xj, ... , x^) = / (xi, X2, . . . ,xif )| 

+ Ph |(xf, xf , ... , x^) 7^ (xi, X2, . . . , Xftr) | /(xf , xf , ... , xj) 7^ / (xi, X2, . . . ,xk)^Ph |/(xf , xJ, ... , x^) 7^ / (xi, X2, . . . , xk)} 



(13) 



Ph |(xf , x°, . . . , x2) ^ (xi, X2, . . . , Xk) \ /(xf , xf, .... x^) = / (xi, X2, . . . , x/f )| 

+ ^ Pif {(xf , X°, . . . , X^) ^ (xi,X2, . . . , Xk) I /(xf ,xj, ... ,X^) = X^} Pif {/(xf , X J , . . . , xj ) = X^ j 



(14) 

P{E} < P |(xf , x°, . . . , x2) 7^ (xi, X2, . . . ,xk) I /(xf , xj, . . . , xJ) = / (xi, X2, . . . , xjc)| 

+ 5Z P {(*f .5?. ■ • ■ ) 5^ (a;i,X2, . . . ,xk) I /(xf ,xj, . . . , x^) = x^} P {/(xf ,xf , . . . , x^^) = x^} . (15) 



a:R7^/(a:i ,X2 ^k) 



P |(xf, xf , ... , Xk) 7^ (Xl, X2, . . . , xk) I /(xf , xf , ... , Xk) = / (Xl, X2, . . . , xk)} 

= P {{toi(xi, X2, . . . , xk) > mi(x'i, X2, . . . , Xk), (x'l, X2, . . . , x'j^) 7^ (xi, ,1,2. . . . , xk)} 

U {mi(xi, X2, . . . , xk) > log{SNR) + m2{x[, x'2, ■ . ■ , x'^), (x'l, x'2, ■ ■ ■ , Xk) 7^ (xi, X2, . . . , xk)} | /(xf , xf , ... , Xk) = / (xi , X2, . . . , xk)} 

(16) 

P |mi(xi, X2, . . . , xk) > mi(xi, X2, . . . , Xk) | /(xf , xf , Xk) = / (xi, X2, . . . , xk)} 



^(xi,x^^,...,x^)es 

(a;^,a:2,... ,a:^ ,a;2 > ■ ■ ■ ,^k) 



+ ^ P |toi(xi, X2, . . . , xk) > log(SJVfl) + m2(x'i, Xj, . . . , Xk) I /(xf , xf , . . . , Xk) = / (xi, X2, . . . , xk)} . (17) 



m4(xi, X2, . . . , XK, Xh) ■ 



'E.ai Xi 



- 2^ hs^oV^bi Xi — HrdV^xr 



+ \oeiSNR). 



(18) 



P |toi(xi, X2, . . . , xk) > m2(xi, Xj, . . . , Xk) + log(SJViJ) | /(xf , xf , . . . , Xk) = / (xi, X2, . . . , xk)} 

= P < mi(xi, X2, . . . , xk) > min mi{x'^ , Xg, . . . , x^, x^) | /(xf , xf , ... , x^) = / (xi, X2, . . . , xk) ^ 

(19) 

< ^ P |mi(xi, X2, . . . , Xk) > m4(xi, X2, . . . , x'^, x^) | /(xf , xf , . . . , x^^) = / (xi, X2, . . . , xk)} . 
«R7^/(«i.»2'--- ■'"k^ 

(20) 

Ph |mi(xi, X2, . . . , xk) > m4(xi , Xj, . . . , x'^, x^) | /(xf , xf , ... , Xk) = /(xi , X2, . . . , xk)} 

K K 

= P/f {I Z£,J^ + IzD^I^ > logiSNR) + \zDi +^/ie,D\/Bra.Ax,|^ + \ZD2 +J2''-s.D^E^biAxi + /ihd v^AxhI^}. (21) 
P /or / * X^-losiSNR) I r - log jSNR) ||x||\ 



logjSNR) , jjxjl 



Tht 

^ = \/Es [hs^D hs2D 



The vector x can be written as 



has a diversity order two. 
ten as. 



as a function of r, attains the minimum value when r = 
^°\iSNH^ ^'^'^ minimum value equals a/2 log{SNR). 
Since, Q{x) is a decreasing function of x, we have. 



aiAxi a2Ax2 
biAxi 62^x2 



ukAxk 
bpcAxK Axfi 



AX 

Since AXAX* is Hermitian, it is unitarily diagonalizable, 
i.e, AXAX* USU*, where U is unitary and S is a 
diagonal matrix. Since AX has a maximum rank two, the 
number of non-zero diagonal entries of S has to be less than 
or equal to two. Let Ai and A2 denote the two diagonal entries 
of S which are possibly non-zero, with Ai > A2. We have 
||x||^ =^SNR hUSU*h*. Let h = hU = [hi h h]- The 
vector h has the same distribution as that of h, since U is 
unitary. 

Since the rank of AX is at least one, Ai > 0. We consider 
the two cases where A2 > and A2 ~ 0. 
Case 1: A2 > 



For this case, upper bounding Q 



logiSNR) 
V2j|x|| 



Q 
Q 



_V2 

log (SNR) 



V2 

we have 



by 



which is upper bounded by e 

< ^-^{\iSNR\hi\''+\2SNR\h2\^) ^23) 



Taking expectation w.r.t \hi\ and \h2\, from ( |23] l, we get, 



< 



(1+ 



Xl SNR 



)(n 



has a diversity order two. 



^ logjSNR) ||x|| 

I. ^/2||x|| + V2 

Case 2: A2 = 

For this case ||x|| — ^/XiSNR\hi\. Hence, 
- \og{SNR) ^ ||x|r 
. 72||x|| V2. 

logjSNR) ^ VXiSNR\hi\ 
V2VMSNR\hi\ V2 



-)■ 



(24) 

Let r = \hi\. Taking expectation w.r.t r, from ( |24| i, we get, 
E (q 



£ 



logjSNR) ||x|| 
V2 



V2\ 



log {SNR) VXiSNRr 



V2VXiSNRr 



^/2 



2re ^ dr 



/ loE(SWH) 
Y X-^SNR 



2Q 



log {SNR) V^iSNRr 



V^VXiSNRr 



+ 



V2 



^ dr 



+ 



L 



XiSNR' 



-2Q 



ii 

log {SNR) 
V2V^iSNRr 



+ 



VXiSNRr 



^/2 



dr . 



In the following, we show that the integrals Ii and I2 



have diversity order two. Note that 



log(SjV-R.) 
V2V)^iSNRr 



V\iSNR.r 
V2 



Q log(SAffl) 
^ IV^^XiSNRt 

we have. 



h < 2Q 



< 



^/\lSNR.r 
V2 



< Q 



^2log{SNR) 



Hence, 



^2\ogiSNR) 



/ log(SJVR) 
V XiSNR 



re ^ dr 



r=0 



SNR 



log(S]VH) 



Since for small x. 



can be approximated as 1 — 



X, at high SNR, we have Ii < 



log (^^]Msm\ 



log(5jV-R) 

SNR \iSNR ■ 



Since 



lim 

SNR^oo log{SNR) 
least two. 



SNR'^ J 



Let To 



' logjSNR) 
\iSNR 

bounded as, I2 < f°° Q 



Let r' 



logjSNR.) 



= 2, Ii has a diversity order at 
The integral I2 can be upper 

log(5jV-R) , VAi5jV-Rr ] , 
^/2^ATS^V7^r ^ V2 J 



^>^^SNR,r ^ ^ function of r, r' 

V 2 



V2s/XiSNRr 

is monotonically increasing for r > tq. Also , for r > rp, r 

^ , r' + x/r'^-2logjSNR) 

can be written m terms or r as, r = ^ 



We have, dr = dr' 

r < 



\/2A7S]VS 



1 



^2XiSNR y ^r'-^-2logjSNR) 

^2\lsNR ' ^2 can be upper bounded in terms of r' as. 



Since 



Hence h < 



Q{r')r' 
^2 log(SiVfl) XiSNR 
1 



1 + 



2log{SNR) 



dr' 



XlSNR J ^2 los(SNR) 

1 21 
1 



Q{r')r'dr' 



Q{r')r'^ 



XiSNR J ^2iog{s NR) \fr^ 



2 \og{SNR) 



-dr . 



Upper bounding Q(r') by e 2 ^ l^i can be shown to be 
upper bounded as ^ snk^ ' which falls as SNR^^. Upper 
bounding Q{r') by e^T", and using the transformation t = 
r'^ — 2\og{SNR), I22 can be upper bounded as, 

1 te-i 



/22 < 



XiSNR^ Jo 



< 



1 



ty/t + 2log{S NR) 
2\og{SNR) 
XiSNR^ io 

1 

e ■ 



+ 2 logl^iV-R) 



XiSNR^ Jo XiSNR^ JO 

2 + 2V2^1og(S'iVi?) 
XiSNR^ ■ 



Vt 



dt 



where the second inequality above follows from the facts 

that ^ = < ^ and ^ = < 1 for 

yt+2iog(sjvfl) - v* ^t+2iog(sjvi?) - 

sufficiently large SNR. The last equality follows from 



the fact that ^dt = V^r(l/2) 
r{z) is the integral, r(z) = 



= V2tt, where 
'H^-^dt. Since, 



lim 



_ , / 2+2v/2¥log(SjVfl) \ 



log(S'7Vi?) 

order 2. This completes the proof of Lemma 1 



2, I22 has a diversity 



Lemma 2: When the 
statement of Theorem 



ability 



two conditions in the 
1 are satisfied, the prob- 
f , . . . ^ {xi,X2, ■ ■ ■ ,xk) I 

^ f{xi,X2,.-. ,Xk), 



The vector x can be written as, 

X = \/Es [hs^D hs^D ■ ■ ■ ^Sa'-D^-R-d] 



aiAxi 02^X2 

^lAcCi 1)2^X2 



uk^xk 
bj^Axx Axfj 



1-2 , . . . , x^) — x'j^} , x'ji 
has a diversity order one. 

Proof: Let 7714 denote the metric as defined in 
( |25| l, given at the top of the next page. Under the 
condition that R transmitted the wrong network coded 
symbol x'^, a decoding error occurs at D only when 

mi{xi,X2,- ■ . ,Xk,x'j^) > m4{x'(,X2,. ■ . ,x'^,x'l^) 
or m4{xi,X2,... ,XK,x'j{) > mi{x'{,X2, ■ . ■ ,x'^), 
for some {x'( , x'^ , . . . , x'l^) ^ {xi,X2, ■ ■ ■ ,xk) 
and x'l^ 7^ /(x'/j X2, . . . , x^). Hence, the 
probability P { (if , xf , . . 



AX 

Since AXAX* is Hermitian, it is unitarily diagonalizable, i.e, 
AX AX* = USU*, where U is unitary and S is diagonal 
with Ai and A2 denoting the two possible non-zero diagonal 
entries, where Ai > Ao- Since the rank of AX is at least one, 
Ai > 0. We have ||x|p = SNR hUSU*h*. Let h = hU = 
[hi h2 hs]. The vector h has the same distribution as that of 
h, since tJ is unitary. Hence, we have, 

P {llxf < loeiSNR)} = P {a,|/.,|^ + X,\h,f < ^°^|^} 



<PiXi\hi\'< 



2 ^ logjSNR) \ 
SNR J ■ 



/(if 



R .t,R 



) 7^ (xi,X2, ■ ■ ■ ,xk) Since, is exponentially distributed. 



c'^} can be upper bounded as in 



( 26) 1 (eqns. (|26| - p9| ) are shown at the top of the next page). 
Using the union bound, from (|26]l, we get (|27]l. 



p[\\k\\^ <\og{SNR)} < ( 



_ log(SWH) 
\ g XiSNR 



Since the matrix 



has rank at 



{x'l 



aiAa;i a2Ax2 ■ ■ ■ ukAxk 
biAxi b2Ax2 ■ ■ ■ BkAxk Axr 
least one for {xi,X2,--- ,xk) 



loejSNR) 

At high SNR, 1 — e >-isnr ^.^jj approximated as 
'Ti'swfl ^ . i^{||x||^ <log(5iVi?)| has a diversity order at 



, 0-2 , • ■ 



where 



Axr 



f{x'{,X2,.-. ,x'^), and Axi — Xi — x'l, the probability 
P{TO4(xi,a;2, . . . ,xk,x'j^) > m4{x'( , X2 , ■ ■ ■ ,x'l^,x'^) 

x^) = x'^} has a diversity order at least 



least one since lim 

SNR^oo 



log 



/ \og[SNR) \ 
\ XiSNR J 



1. Since Q{x) < 



I f{xi,x 
one. 



R 

2 : ■ 



Ph {m4{xi,X2, . ■ ■ ,xk,x'j{) > mi{x'l,x'^,... ,x'^) 



R ^R 
^2 1 ■ 



. . ,x^) = x'^} can be written in terms of the 



I /(if 

additive noise and ZD2, as given in (|28 

Let xi = Y.f=i hs.DCiiy/EsXt, X2 = '^'^ 
hRDAxR)y/El. Also, let zd = [zot 



log {SNR) 

e~T" ^ the integral on the right hand side of ( |31[ ) can be upper 
bounded as, 

\og{SNR) ^ ||x|r 
V2. 



(I]i=i hsiohxi 



H:{\\^f>\og(SNR)} 

<[ 

J H:{\\xf>\og{SNR)} 



V2||x|| 

l~ llxll - 



dH 



-dH. 



(32) 



ZD2 



and 



[a:;i X2] ■ Then ( [28] l can be simplified as in ( |29| l, where 
w = V2i?e{zQ^}, is distributed according to A/'(0,1). 
Hence, 



log (S'A^i?) 



1 



{\\xr<log{SNR)} 



V2 

1-Q 



We consider the following two cases when A2 > and A2 
0. 

Case 1: A2 > 

For this case, from the integral in ( (32] i, we get, 

log(5'iV_R) II X 



H:{||x||2>log{S]V_R)} 



^/2| 



^/2 



dH 



+ l{||xf >log(SJVfl)}Q 



log (SNR) 
V2\\k\\ 
log (S'A^i?) 
V2\U\\ 



l|xH 



+ 
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(30) 



lhlP=0-'|h2p=0 



hir„-|h2|- 



d\hi\^ d\h2\ 



Taking expectation with respect to the fade coefficients in 



30)1, P{ 



^ ^ log(SjV-R) _ Jlxjl ]. 
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log (SNR) 
^llxll 



^ I can be upper bounded as. 



(l + ^if^)(l+^^)' 

which falls as SNR-^ at high SNR. 
Case 2: A2 = 



< 



H:{||xl|">log(SAfi?.)} 



P{||x||' <log(5Ari?)} 
log {SNR) 



V2 



dH. 



(31) 



m4,(xi ,X-2. . . . .Xk, Xft) ^ 



- \og{SNR). 



(25) 



P {{xi ,x" , . . . ,x'^) =t (xi,X2., . . . ,xk) \ /(if, if , . . . , = 

= P {{vai{xi,X2, . . . ,xk,x'j^)> mi{x'J^, x'g , x'^) , x'^ ^ /(^^a, ^b)} U {m4,{xi, X2, . . . ,xk,x'„) > mi (k^ , Kg ) } | /(if , if , . . . , if ) = a^j 



(x^,x^);i(xi,X2,... ,xj^) 



< 



P {m4(a;i, 2:2, ■ • ■ , kk, a:^,) > "14(2;^, a;^ , x'^) | /(if , if , . . . , if ) = x^j 

(x^,x^)es2 x'^es. 

(x^,x^)^(xi,X2,.^. ,x^) x'^^/(x'^,x'^) 

+ -P{m4(2:i, 2:2, ■ ■ ■ , , x'^) > mi(2;'^, | /(if 

^(x^,x'^)es2 

(X^,X^)^(X4,X2,... :X^) 



(26) 



^k) = ^'n] ■ (27) 



Pff |log(SWP) + m4(xi, 2:2, ■ ■ ■ ,XK,x'fi) > mi(x'J^,x'g) \ /(if , if , . . . , if ) = 2;^ | 

K K 

= Ph{\ zo^f + \zD2f +log(SWP) > \zDi +^hs,DVE^aiAx, + \zd2 +J2'^SiD^/lhb,Ax, + hnoAxR\^}. (28) 



log (SNB) 



log{SNR) ||x|| 
y2||x|| V2 



(29) 



For this case. 



if:{||xlp>log(SArp)} 



f{H)dH 



< 



\hi\> 



■y-lasjSNR) 
y/^l SNR 



Q 



logjSNR) ^ M 
V2||x|| V2 
VXiSNR\hi\ - y/\og{SNR) 
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2|/ii|e-l''il (33) 



The above inequality follows from the fact that for \hi\ > 



Let r = From ( [33] l, since (5(a;) < e 2 the integral / 
can be upper bounded as. 



/ < 



XiSNR\ 



■yXiSNRlogiSNR) 
4 + \iSNR 



VAi SNR 



fc2 



Since J e-'=i('^-'=^)'rfr = Q [(rp - fc2)V2fci 



/ < J^g-fci(ro-fc2)^ 
- 2fci 



^2 k3\l -j^Q 



Since Q [(ro - fc2)y2l];] < 1, e'''^^'''''''^^^ < 1, and fcg 
can be approximated as one at high SNR, substituting for 
To, fci, k2 and k^, we get, 

1 v/'^AiS'A^i?log(57Vi?) 



/ < 



2(l + ^if^) 



(4 + Ai57Vi?) Jl + 



Since the above upper bound on / falls as SNR^^ at high 
SNR, I has a diversity order at least 1. This completes the 
proof of Lemma 2. ■ 



log(SiV«) I log(SWfl)AiSJVH 
g 4 ^ 16+4AiSJVfi (^j. 



fe3 



(34) 



Let ro, fci, k2 and fca be defined as shown in ( |34| l. From 
(|34|i, the upper bound on / can be written as. 



I<—[ 2(r - fc2)e-'=i('-'=^''dr + k2 ks f e''"'^''-''^''^ dr. 



